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ABSTRACT 

The TAPPS employs the MASS model to produce 
mesoscale atmospheric simulations in support of the 
Wake Vortex project at Dallas Fort-Worth 
International Airport (DFW). A post-processing scheme 
uses the simulated three-dimensional atmospheric 
characteristics in the planetary boundary layer (PBL) 
to calculate the turbulence quantities most important to 
the dissipation of vortices: turbulent kinetic energy and 
eddy dissipation rate. TAPPS will ultimately be 
employed to enhance terminal area productivity by 
providing weather forecasts for the Aircraft Vortex 
Spacing System (AVOSS). The post-processing scheme 
utilizes experimental data and similarity theory to 
determine the turbulence quantities from the simulated 
horizontal wind field and stability characteristics of the 
atmosphere. Characteristic PBL quantities important to 
these calculations are determined based on 
formulations from the Blackadar PBL 
parameterization, which is regularly employed in the 
MASS model to account for PBL processes in 
mesoscale simulations. The TAPPS forecasts are 
verified against high-resolution observations of the 
horizontal winds at DFW. Statistical assessments of the 
error in the wind forecasts suggest that TAPPS 
captures the essential features of the horizontal winds 
with considerable skill. Additionally, the turbulence 
quantities produced by the post-processor are shown to 
compare favorably with corresponding tower 
observations. 
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INTRODUCTION 

The Terminal Area PBL Prediction System 
(TAPPS) is being developed to produce meso-p scale 
simulations of meteorological conditions in support of 
the Wake Vortex project at Dallas Fort-Worth 
International Airport (DFW ). 1 The centerpiece of 
TAPPS is the Mesoscale Atmospheric Simulation 
System (MASS) 2 , a mesoscale numerical weather 
prediction model that provides real time, operational 
simulations of the planetary boundary layer (PBL) 
conditions in the DFW area. TAPPS will ultimately be 
employed to enhance terminal area productivity by 
providing weather forecasts for the Aircraft Vortex 
Spacing System (AVOSS ). 3 ' 4 The Wake Vortex project 
in general and AVOSS in particular is designed to, at a 
given terminal, assess the wake vortex drift and 
dissipation under varying weather conditions, and 
compute the optimal safe spacing for landing 
commercial aircraft. 

This study will focus upon the overall performance 
of the TAPPS system as well as presenting recent 
enhancements in the post-processing system. The post- 
processor has been developed by the Department of 
Marine, Earth, and Atmospheric Science at the North 
Carolina State University to compute the predictive 
quantities from the TAPPS simulations required by 
AVOSS. The scheme employs the simulated three- 
dimensional winds, thermal structure, and moisture 
characteristics of the atmosphere to calculate the 
turbulence-related variables considered most important 
to the dissipation of vortices: the turbulent kinetic 
energy (TKE) and eddy dissipation rate (EDR). 




MASS/TAPPS MODEL CHARACTERISTICS 


Hydrostatic primitive equation model 
3-D primitive equations for u, v, T, q, and p 
Cartesian grid on a polar stereographic map 
Op terrain-following vertical coordinate system 
One-way interactive nested grid system 
Time-dependent lateral boundary conditions 
First guess provided by large scale gridded analyses 
Reanalysis of rawinsonde, surface, and asynoptic data 
using a 3-D optimum interpolation scheme 
Blackadar PBL scheme 
Surface energy budget 

Prognostic equations for cloud water, ice, rain water, 
and snow 

Convective parameterization scheme 
Post-processing system to compute AVOSS input 
quantities 


POST-PROCESSING SYSTEM 


In order for AVOSS to effectively predict the drift 
and dissipation of vortices in the approach corridor of 
DFW, it is necessary for the post-processor to provide 
highly-detailed forecasts of winds, particularly the 
cross-runway (U) component, from the TAPPS 
simulations. Time-height cross-sections of the winds 
are computed by the system with a temporal resolution 
of 30 minutes, and a vertical spacing of 15m between 
the surface and 1000m. The variance and the shear of 
these winds are also very important for AVOSS, and 
must be computed with precision. Finally, in order to 
assess vortex dissipation, diagnostics are computed that 
evaluate the growth, dissipation, and magnitude of the 
turbulence. 

The current version of the TAPPS post-processor 
employs the model winds, temperature, and moisture 
fields in conjunction with a first-order closure scheme 
developed by Han et al. to calculate the EDR and TKE 
for AVOSS. 1 The scheme utilizes a series of equations 
developed from experimental data and similarity 
scaling to determine the EDR and TKE from 
observations of the three-dimensional wind field and 
the moist thermodynamic stability characteristics of the 
atmosphere. 5 The scheme assumes that the PBL has two 
regimes: 1) neutral and stable, and 2) unstable. 

For a statically stable or neutral PBL, the scheme 
divides the results into two layers. For the surface layer, 
which is defined as ( Z < 0.1 h ), where Z, is the height 
above the ground and h is the PBL height, the EDR 
( £ ) is given by: 
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where M* is the friction velocity, k ~ 0.4 is the von 

Karman constant, and L is the Monin-Obukhov length 
scale. Above the surface layer. 
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For a statically unstable PBL, in which a convective 
regime dominated by buoyancy prevails, the EDR is 
calculated as follows: 
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for the surface layer ( z < 0.1 h). Above the surface 
layer in a convective regime, turbulent processes arising 
due to the statically unstable conditions lead to well- 
mixed atmospheric conditions up through the top of the 
PBL. In this so-called mixed layer , the EDR becomes: 
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where W s , is the convective velocity scale. 

The formulation for the TKE (e) follows a similar 
framework. For a statically stable or neutral PBL, the 
scheme divides the results into two layers. For the 
surface layer ( z ^ 0. lh ), 


e = 6 ul , 
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and above the surface layer. 
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For a statically unstable PBL, the TKE is calculated as 
follows: 
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for the surface layer ( 7 < (). 1/7 ), and: 


regime, (3) the forced convection regime, and (4) the 
free convection regime. These regimes are determined 
by calculating the bulk Richardson number and 
assessing the buoyant characteristics of the PBL. Then 
U t becomes: 
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for the mixed layer. 

There are three quantities in the above equations 
that must be defined before this scheme can be applied 
to the simulated fields: the PBL height ( h ), the friction 
velocity ( U t . ), and the convective velocity scale ( W* ). 
In previous versions of the post-processor, these 
quantities were derived by first evaluating the turbulent 
character of the PBL from the resolved model winds 
and then inferring the other quantities based on as series 
of empirical relationships between the TKE and the 
generally observed characteristics of the PBL. 
However, Eqns. 1-8 were designed to calculate the 
EDR using wind values taken from tower observations 
of PBL winds. The resolved winds from a mesoscale 
atmospheric model (which has a horizontal resolution 
on the order of 12km) will not, in general, display the 
same degree of small-scale variability that occurs in 
tower-based PBL observations. As a result, the model 
resolved winds will not provide information about the 
turbulent characteristics of the PBL, and thus will tend 
to yield values for EDR and TKE that are significantly 
smaller than those calculated from observed winds. 

In order to overcome this potentially serious 
limitation, the post-processing scheme has been revised 
such that the EDR and TKE calculations are more 
representative of the background turbulence than the 
resolved mesoscale model fields. This was 
accomplished by using formulations from the Blackadar 
PBL parameterization, which is employed within 
MASS to account for sub-grid scale PBL processes that 
cannot be explicitly calculated by the model due to its 
resolution. 6 7 By using the Blackadar parameterization, 
representations for W* , h , and tv, that are not wholly 
dependent upon the resolved model winds can be 
obtained. 

Perhaps the most important of these modifications 
to the overall performance of the post-processor is the 
formulation for W* . Within the Blackadar PBL 
parameterization, W* is calculated by using a scheme 
that separates PBL conditions into 4 regimes: (1) the 
stable regime, (2) the damped mechanical turbulence 


V 


where V a is the velocity at the top of the surface layer, 
Z a is the height of the surface layer, Z a is the roughness 
length, and 'F is the correction to the log profile of 
momentum due to deviations from neutral stratification 
within the various PBL regimes. W* is a pre -defined 

minimum of the friction velocity that is generally used 
only in the stable regime (regime 1). This 

parameterization for U t . , which takes into account the 
background thermodynamic and turbulence 
characteristics of the simulated atmosphere, provides 
substantially better information about the turbulent 
characteristics of the PBL, which in turn leads to more 
accurate and physically consistent EDR and TKE 
predictions. 

The Blackadar PBL parameterization also includes 
a formulation for the PBL height ( h) in the convective 
regime (regime 4) based on the buoyancy 

characteristics of the airmass. 6 ' 7 Outside of the 
convective regime, an empirical relationship that 
defines h for a stationary, neutral PBL is employed. 5 

The convective velocity scale is represented by: 
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where g is the gravitational acceleration, T is the 

reference temperature, and {wd ) is the mean surface 

heat flux. Since the surface heating parameterization 
within the MASS model accounts for the shortwave 
radiative heating effects, a representation of the surface 
heat flux that is physically consistent with the W* and 

h calculations can be used. By incorporating this value 
into Eqn. 10, the convective velocity scale will likewise 


3 

American Institute of Aeronautics and Astronautics 



be consistent with the rest of the input into the EDR and year presents a considerable challenge for AVOSS, and 
TKE scheme. that TAPPS regularly captures these features is also 


encouraging. Furthermore, that TAPPS is providing 



TAPPS RESULTS 


In order to assess the performance of the post- 
processor, it is necessary to first verify that the model 
winds are in agreement with observations. Figs. 1-2 
show comparisons between the TAPPS simulated U 
(cross-runway) and V (along-runway) components of 
the horizontal winds and the wind consensus profiles 
produced by MIT Lincoln Labs from tower, profiler, 
lidar, and other in situ observations at DFW (referred to 
as the AVOSS Winds Analysis System, or AW AS). 8 
Fig. 1 shows time-height cross sections of the two wind 
components for 2 November, 1999, while Fig 2 
corresponds to 19 November, 1999. These two cases 
were chosen because both situations involve the 
passage of a synoptic-scale front, which is traditionally 
a difficult feature for mesoscale models to predict. As 
evidenced by the sharp reversal in direction of the V- 
wind components (Figs. lc,d and 2c, d), TAPPS and the 
observations both clearly indicate the passage of the 
fronts in these cases. 

It is noteworthy that TAPPS reproduces with some 
degree of accuracy both the level and magnitude of the 
wind maximums associated with the approaching 
fronts. In particular, the agreement between TAPPS and 
AW AS for the U-wind component maximum on 2 
November (Fig. la,b) and for the V-component 
maximum on 19 November (Fig. 2c, d) shows 
considerable skill. The high degree of agreement 
between the overall structures in both wind components 
during the passage of these fronts suggests that TAPPS 
can reproduce mesoscale features in the horizontal wind 
fields that are often difficult to replicate. Although not 
shown here, it is also notable that TAPPS also verifies 
well in situations where a nocturnal low-level jet (LLJ) 
is observed by the AW AS. The almost ubiquitous 
nature of LLJs in the DFW area during much of the 


A statistical error analysis that assesses the overall 
performance of TAPPS in comparison to AW AS is 
shown in Fig. 3. These statistics were computed from 
the differences between TAPPS and AWAS for 28 real- 
time simulations in November 1999. These simulations 
were initiated at 0000 UTC on each date. Results were 
analyzed for every 30-minute time period for which 
there were both observations and model output, 
resulting in a sample size on the order of 1300 for each 
level. For the purpose of understanding these results, it 
is worth noting that the estimated error in the 
instruments that provided data for the AWAS profiles is 
on the order of 1 m/s. 8 The bias statistics suggest that 
there is little overall bias in the TAPPS predictions, 
particularly in the lowest levels. Even at the highest 
levels, the overall bias is less than the estimated errors 
in the observations. 

The mean absolute error (MAB) is on the order of 
1.5 m/s in the lowest 500m of the profile, and increases 
to near 2 m/s at the higher levels. The increase in MAB 
at the higher levels can be attributed to the typically 
higher wind speeds in the upper levels, which causes 
larger magnitude differences should the model develop 
phase errors in the simulated weather systems. 
Nevertheless, overall MAB values on the order of 1.5 
m/s in the lowest 500m, which is the layer most vital to 
the performance of AVOSS, is very encouraging. The 
RMS errors are, as expected, somewhat larger than the 
mean absolute errors. However, since RMS error 
calculations are designed to punish particularly large 
differences between TAPPS and AWAS, RMS errors 
less than 1 m/s greater than the MAB in all cases 
suggests that there are few extremely large "misses" 
during the month. 

Finally, Figs. 4 and 5 show an example of TAPPS 
TKE and EDR predictions respectively at three 
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different times on 15 November, 1999. The TAPPS 
profiles are compared against EDR and TKE profiles 
inferred from 5m and 45m tower observations by using 
the technique described in Han et al. 5 These profiles 
show that, in this particular case, TAPPS generated 
TKE and EDR values that are of the same order of 
magnitude as those inferred from the tower 
observations, which is very encouraging. The vertical 
structure of the quantities is also well represented in the 
early times in particular (Figs 4a,b and 5a,b), when the 
PBL was in a stable regime. In the afternoon (4c and 
5c), the vertical structure of the EDR and TKE are less 
comparable, due to disagreements between the TAPPS 
simulations and the tower observations of EDR and 
TKE in the convective regime. Nevertheless, the fact 
that the two techniques agree within an order of 
magnitude shows a considerable improvement over the 
techniques that were previously employed as part of the 
TAPPS system. However, there is clearly a need for 
more comprehensive statistical comparisons between 
TAPPS and observations of EDR and TKE. There are 
plans to start calculating the EDR and TKE by using 
AW AS and profiler data in the near future, which will 
provide more comprehensive verification information 
with which the overall performance of the post- 
processor can be assessed. 

CONCLUSIONS 

A new formulation for predicting the evolution of 
turbulence variables in the PBL from TAPPS 
simulations has been presented. The physics of the new 
formulation, which utilizes elements of the Blackadar 
PBL parameterization to account for the sub-grid scale 
characteristics necessary to reproduce the details of 
turbulent flow, was described. 

The overall performance of TAPPS was assessed 
by comparing TAPPS and AW AS horizontal winds in 
the lowest 1000m of the atmosphere. Examinations of 
two specific cases as well as an overall statistical 
assessment of the differences suggested that TAPPS 
reproduces the salient PBL characteristics in the 
horizontal wind fields with considerable skill. 
Comparisons between TAPPS simulation results from 
the redesigned post-processor and tower-based 
observations of EDR and TKE showed that the new 
formulation agreed well with observations for one 
particular case. Both the temporal evolution and the 
magnitudes of the simulated quantities were shown to 
be representative of the conditions at DFW. 
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TAPPS 



AWAS 


Figure 1: Time-height cross section of the horizontal wind components for 2 November, 1999 at DFW in 
m/s. a) TAPPS U-component. b) AWAS U-component. c) TAPPS V-component. d) AWAS 
V-component. 
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Figure 1: Continued 
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Figure 2: Time-height cross section of the horizontal wind components for 19 November, 1999 at DFW in 
m/s. a) TAPPS U-component. b) AWAS U-component. c) TAPPS V-component. d) AWAS 
V-component. 
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Figure 3: Vertical profiles of bias, mean absolute error, and rms errors between TAPPS and AWAS wind 
components for the month of November in m/s. a) U-wind component, b) V-wind component. 
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c) 


1900 UTC 


TAPPS vs. Tower-derived TKE 



Figure 4: Continued 
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1900 UTC 


TAPPS vs. Tower-derived EDR 



Figure 5: Continued 
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